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Abstract 

Neutrino oscillations are studied in the general framework of open quan- 
tum systems by means of extended dynamics that take into account pos- 
sible dissipative effects. These new phenomena induce modifications in 
the neutrino oscillation pattern that in general can be parametrized by 
means of six phenomenological constants. Although very small, strin- 
gent bounds on these parameters are likely to be given by future planned 
neutrino experiments. 
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1. INTRODUCTION 



A large variety of open quantum systems can be modeled as being subsystems in in- 
teraction with large environments. The time evolution of the total system, subsystem plus 
environment, is unitary and follows the standard rules of quantum mechanics. However, 
the dynamics of the subsystem alone, obtained by eliminating the environment degrees of 
freedom, is no longer unitary, as it develops dissipation and irreversibility. [1-3] 

When there are no initial correlations between subsystem and environment and their 
mutual interactions can be considered weak, the resulting subdynamics can be described 
in terms of so-called quantum dynamical semigroups. These are time-evolution maps that 
encode very general physical requirements, like entropy increase, forward in time compo- 
sition law (semigroup property) and complete positivity; these properties are essential for 
the correct physical interpretation of the subdynamics. 

Although this description of open quantum systems has been originally developed 
in the framework of quantum optics, [4-6] it is very general and can be applied to model 
a variety of different phenomena. Recently, it has been adopted to study the effects of 
dissipation and irreversibility in various particle physics phenomena. [7-14] 

The original motivations for such investigations were based on quantum gravity effects: 
[15-20] due to the quantum fluctuations of the gravitational field and the appearance of 
virtual black holes, spacetime becomes "foamy" at Planck's scale, leading to possible loss of 
quantum coherence. From a more fundamental point of view, also string theory could lead 
to similar effects in the low energy domain. [21] Indeed, subdynamics described by quantum 
dynamical semigroups are the result of the interaction with a "gas" of DO-branes at Planck's 
temperature, obeying infinite statistics. [7] Nevertheless, since not enough details about the 
"microscopic" dynamics are known to allow precise estimates of the magnitude of these 
new effects, the description of dissipative phenomena that we shall discuss below should 
be thought as being phenomenological in nature. 

These new, non-standard effects are very small, since they are suppressed by inverse 
powers of the Planck mass, as a rough dimensional estimate suggests, and therefore very 
difficult to observe in practice. However, there are particular situations, involving inter- 
ference phenomena, in which they might be in the reach of present or future experiments. 
Indeed, detailed studies of neutral meson systems, [8, 10, 13] and neutron interferome- 
try,[14] using quantum dynamical semigroups have already been performed and order of 
magnitude limits on some of these dissipative effects have been derived using available 
experimental data. [9, 11, 14] One of the most interesting outcome of these investigations 
is that future experiments, in particular those involving correlated neutral mesons, should 
be able to ascertain with high accuracy the presence of such dissipative phenomena. 

Neutrino physics is certainly another obvious place where to look for non-standard 
effects. Many neutrino experiments are presently taking data and other will start operating 
in the near future, so that it appears timely to discuss in detail to what extent dissipation 
can affect those observations. 

We shall limit our considerations to the vacuum oscillations of two species of neutrinos. 
In this case, possible dissipative effects can be parametrized in terms of six phenomeno- 
logical constants that modify the pattern of the transition probability V among the two 
neutrino flavours, by introducing exponential dumping factors. Although the explicit ex- 
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pression of V is in general rather complicated, in the generic case its asymptotic (large 
time) behaviour turns out to be independent from the mixing angle. Various approximated 
expressions for V will also be discussed; they can be of help in fitting the experimental data. 
Finally, in the last section we shall present a discussion on a possible physical mechanism 
that could be at the origin of the dissipative phenomena. 



2. QUANTUM DYNAMICAL SEMIGROUPS AND NEUTRINO 
OSCILLATIONS 

Quite in general, states of a quantum system evolving in time can be described by a 
density matrix p; this is a hermitian, positive operator, i.e. with positive eigenvalues, and 
constant trace. We shall analyze the evolution of neutrinos created in a given flavour by 
the weak interactions and subsequently detected at a later time. Assuming the neutrinos to 
be ultrarelativistic, the study of the transition probability for the original tagged neutrinos 
to be found in a different flavour can be performed using an effective description; [22-24] 
further, for simplicity, we shall limit our considerations to the mixing of two neutrino 
species. Then the neutrino system can be modeled by means of a two-dimensional Hilbert 
space, taking as basis states the two mass eigenstates. 

With respect to this basis, the two flavour states, that conventionally we shall call 
V e " and V M " , are represented by the following 2x2 matrices: 

/ cos 2 9 cos 9 sin 9 \ , , 

Pv * ~ Vcos#sin# sin 2 # J ' {ZAa) 

( sin 2 9 — cos 9 sin 9 \ , . 

p ^"V-cos#sin# cos 2 # y= A -P"e> 

where 9 is the mixing angle. 

As explained in the introductory remarks, our analysis is based on the hypothesis that 
the evolution in time of the neutrino state p is given by a quantum dynamical semigroup, 
i.e. by a completely positive, trace-preserving family of linear maps: p(0) i— > pit). These 
maps are generated by equations of the following form: 

^ = -iH eS p(t) + ip(t) tfeff + L[p{t)\ . (2.2) 

The first two terms in the r.h.s. of this equation are the standard quantum mechanical 
ones, that give rise to the traditional description of neutrino oscillations. They contain the 
effective (time-independent) hamiltonian if e ff; neglecting effects due to possible neutrino 
instability, it can be taken to be hermitian. The additional piece L[p] is a linear map, whose 
form is completely fixed by the conditions of complete positivity and trace conservation: [1] 

L\p] = -\T l {A]A j p + pA]A j ) +$> J M], (2.3) 



3 



where the operators Aj must be such that ^ ■ AjAj is a well-defined 2x2 matrix. The 
additional requirement of entropy increase can be easily implemented by taking the Aj to 
be hermitian. [8] It should be stressed that in absence of L[p], pure states (i.e. states of 
the form \ifj) (if;\) would be transformed into pure states. Only when the extra piece L[p] 
is also present, p(t) becomes less ordered in time due to a mixing-enhancing mechanism: 
it produces dissipation and irreversibility, and possible loss of quantum coherence. 

As already mentioned, equations of the form (2.2), (2.3) have been used to describe 
various phenomena related to open quantum systems; in particular, they have been applied 
to analyze the propagation and decay of neutral meson systems. [8-13] Although the basic 
general idea behind these treatments is that quantum phenomena at Planck's scale produce 
loss of phase-coherence, one should keep in mind that the form (2.2), (2.3) of the evolution 
equations is independent from the microscopic mechanism responsible for the dissipative 
effects. Indeed, it is the result of very basic physical requirements that the complete 
time evolution, j t : p(0) i— > p(t), needs to satisfy; generally, the one parameter (=time) 
family of linear maps 7< should transform density matrices into density matrices and have 
the properties of increasing the von Neumann entropy, S = — Tr[plnp], of obeying the 
semigroup composition law, 7t[p(t')] = p(t + t'), for t, t' > 0, of being completely positive. 
[1-3] In view of this, the equation (2.2), (2.3) can be regarded as phenomenological in 
nature; nevertheless, possible physical mechanisms leading these equations will be discussed 
in the final section. 

Among the just mentioned physical requirements, complete positivity is perhaps the 
less intuitive. Indeed, it has not been enforced in previous analysis, in favor of the more 
obvious simple positivity. Simple positivity is in fact generally enough to guarantee that 
the eigenvalues of the density matrix p(t) remain positive at any time; this requirement is 
obviously crucial for the consistency of the formalism, in view of the interpretation of the 
eigenvalues of pit) as probabilities. 

Complete positivity is a stronger property, in the sense that it assures the positivity of 
the density matrix describing the states of a larger system, obtained by coupling in a trivial 
way the neutrino system with another arbitrary finite-dimensional one. Although trivially 
satisfied by standard quantum mechanical (unitary) time-evolutions, the requirement of 
complete positivity seems at first a mere technical complication. Nevertheless, it turns out 
to be essential in properly treating correlated systems, like two spin-zero neutral mesons 
coming from the decay of a vector-meson resonance; it assures the absence of unphysical 
effects, like the appearance of negative probabilities, that could occur for just simply 
positive dynamics. [25] For these reasons, in analyzing possible non-standard, dissipative 
effects even in simpler, non correlated systems, the phenomenological equations (2.2) and 
(2.3) should always be used.^ 

In the case of the neutrino system, a more explicit description of (2.2), (2.3) can be 
given. In the chosen basis, the effective hamiltonian that gives rise to the standard vacuum 

t We have argued before (see also the discussion in Section 5) that the microscopic 
mechanism leading to the non-standard, dissipative phenomena are likely to originate from 
quantum gravity or string effects. They presumably act in the same way for all systems; 
it is therefore unjustified to adopt different formulations for correlated and uncorrelated 
systems. 
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oscillations can be written as: [22-24] 



*-=rr (2 - 4) 

where i? is the average neutrino energy, while uj = Am 2 /4i? encodes the level splitting due 
to the square mass difference Am 2 of the two mass eigenstates. In the case of oscillations 
in matter, H e e has a more complicated expression, that takes into account the coherent 
interactions of the neutrinos with the matter constituents. For simplicity, in the following 
we shall limit our discussion to vacuum oscillations: we are in fact interested in studying 
possible dissipative effects, which are quite independent from the specific form of the 
standard effective hamiltonian. 

The explicit expression of the term L[p] in (2.3) can be most simply given by expanding 
the 2x2 matrix p in terms of Pauli matrices (7j and the identity a®: p = p^o M , p = 0, 
1, 2, 3. In this way, the map L[p] can be represented by a symmetric 4x4 matrix [L^ u ~\ , 
acting on the column vector with components (po, p\, p2, /O3). It can be parametrized by 
the six real constants a, 6, c, a, (3, and 7: [8] 



[L^v] = -2 



/0 

a b c 

b a (3 

\0 c 7. 



(2.5) 



with a, a and 7 non-negative. These parameters are not all independent; the condition of 
complete positivity of the time-evolution p — > p{t) imposes the following inequalities: 

2R = a + -f-a>0 , U = RS - b 2 > , 

2S = a + 7- a>0, V = RT - c 2 > , 

— — (2 6) 

2T = a + a-7>0, Z = ST - [3 2 > , 

X = RST - 2 bc(3 - Rf3 2 - Sc 2 - Tb 2 > . 

Taking into account that the equation in (2.2) is trace preserving, from the initial 
normalization condition Tr[p(0)] = 1, one immediately obtains that po = 1/2, for all 
times. Then, the evolution equation for the remaining three components of pit) can be 
compactly rewritten in a Schrodinger-like form: 

^|p(f)> = -2W|p(t)> (2.7) 

where the vector \p(t)) has components (pi, P2, P3), and 

/ a b + uj c\ 
H= b-u a p\. (2.8) 
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The formal solution of (2.7) involves the exponentiation of the matrix H: 

\p(t)) = M(t) |p(0)> , M(t)=e- 2nt . (2.9) 

Let us assume that at the beginning the neutrinos were of type V e " ; the probability 
of having a transition into the type "z/ M " at time t is given in our formalism by: 

V„ e ^(t) = Tr[p„ e (t) p^} , (2.10) 

where Pu e (t) is the solution of (2.7) with the initial condition given by the matrix p l/e in 
(2.1). Using (2.9) and the matrices in (2.1), one explicitly finds: 

V„ e ^(t) = 1 1 cos 2 29 [l-M 33 (t)] + sin 2 20 [l - M n (t)] - 1 sin 40 [M 13 (t) + M 31 (t)] j . 

(2.11) 

When the additional piece L[p] in (2.3) is not present, one simply obtains: 

M n (t) = cos(2c^) , M l3 (t) + M 31 (t) = , M 33 (t) = 1 , (2.12) 

so that (2.11) reduces to the well known standard expression for the oscillation probability 
in vacuum: 

Vi%^(t) = sin 2 26 sin 2 ^. (2.13) 

Therefore, any deviation from (2.12) that might be found in fitting the expression (2.11) 
with data from neutrino experiments would provide evidence for dissipative phenomena in 
neutrino physics. ^ 



3. TRANSITION PROBABILITY: GENERAL PROPERTIES 

Explicit expressions for the entries of the matrix M(t) appearing in (2.9) can be 
given by diagonalizing 7i in (2.8); this can always be done by solving the corresponding 
eigenvalue equation, 

H\v^) = A (fc) \v ik) ) , k = 1,2,3, (3.1) 

via Cardano's formula. [26] Then, using the diagonalizing matrix [Dgk] = v^ k \, built with 
the components of the eigenvectors \v^), one formally writes: 

3 

M lJ (t) = J2e- 2X(k)t D lk D~} . (3.2) 
k=i 

t Different physical mechanisms have been proposed in the literature to account for the 
observed neutrino flux deficit: they all predict expressions for the transition probability 
V Vk ^ v (t) that differ from that in (2.11); see the discussion at the end of Section 5. 
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The explicit expressions of A^ fc ^ and [D^] in terms of the dissipative parameters a, 
6, c, a, 7 and uj is however cumbersome, making the formula (3.2) unmanageable 
in practice; for this reason, we shall discuss particularly interesting limits of the general 
expression (2.11) for the transition probability V Vf ,^ v ^{f) in the next section. Nevertheless, 
general conclusions on the behaviour of (3.2) can be obtained by studying in more detail 
the eigenvalue problem in (3.1). 

The three eigenvalues A^, X^ 2 \ of the matrix H are solutions of the cubic equa- 
tion: 

A 3 +rA 2 + sA + w= 0, (3.3) 

with real coefficients, 

r = -(A (1) + A (2) + A (3) ) = -(a + a + 7) , (3.4a) 



s = + A^A^ + A^ 3 ) = aa + a 7 + a 7 " b 2 - c 2 - (3 2 + uj 2 , (3.46) 



w ee -A (1) A (2) A< 3 > = a(3 2 + ac 2 + 7 (o 2 - w 2 ) - aa 7 - 2 bc(3 . (3 Ac) 

According to the sign of the associated discriminant V = p 3 + g 2 , p = s/3 — (r/3) 2 , 
q = (r/3) 3 — rs/6 + w/2, the eigenvalues are either all real (T> < 0), or one is real and the 
remaining two are complex conjugate (T> > 0). The degenerate case T> = occurs when 
two real eigenvalues are equal; all three coincide for p = q = 0. 

Furthermore, the quantum dynamical semigroup generated by (2.2), (2.3) is bounded 
for any £,[27] so that the real parts of A*- 1 -*, A^ 2 -*, A^ are surely non-negative (otherwise 
the entries My(t) in (3.2) would blow up for large times). 

When u = 0, the matrix H is real, symmetric and non-negative, as guaranteed by the 
inequalities (2.6); therefore, its eigenvalues are all real and non-negative: V < and this is 
possible only for p < 0. The discriminant V starts becoming positive only for sufficiently 
large uj, since, as it is clear from the definitions (3.4), the contribution of uj to p is equal 
to cu 2 /3, and thus it is positive. 

Therefore, the time-behaviour of the transition probability V Ue ^ Ull (t) depends on the 
relative magnitude of uj with respect to the non-standard parameters a, 6, c, a, (3 and 7. 
In particular, an oscillatory behaviour is possible only when the dissipative parameters are 
small compared to a;; on the other hand, when dissipation is the dominant phenomenon, the 
time-dependence in (3.2), and therefore in (2.11), is characterized by exponential dumping 
terms. 

This analysis allows a general discussion on the asymptotic behaviour of V Vf ,^ Vlx {$) 
for large £.[27, 1] In the generic case, det(H) 7^ and all three eigenvalues A^, \( 2 \ A^ 
are thus non-vanishing, with positive (or zero) real part, as discussed above. When T> < 0, 
the eigenvalues are all real, so that all entries of the matrix M(t) in (3.2) approach zero 
for large £, due to the exponential dumping factors. The same is true also in presence 
of two complex conjugate eigenvalues, unless their real part is identically zero. However, 
this situation never occurs when there is a non-vanishing dissipative contribution (2.5) in 
the equation (2.2). Indeed, from (3.4) one finds that the condition for having two purely 
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imaginary eigenvalues is given by: w — rs = 0; recalling the definitions in (2.6), it can be 
rewritten as: X + (R + S + T) [U + V + Z] + 2 (R + S + T) 3 + u 2 (R + S + 2T) = 0. Since 
by the inequalities in (2.6) all the terms in the l.h.s. are non negative, they must be zero 
separately, which is possible only for a = b = c = a = (3 = ^= 0. 

Therefore, in presence of dissipative phenomena, the generic large t behaviour of the 
transition probability in (2.11) is independent from the mixing angle 9: 

V^.it) ~ \. (3.5) 

The situation might be different however when det(7i) = and we are in presence of 
zero eigenvalues. In this special case, u> and the dissipative parameters a, 6, c, a, (3 and 7 
need to satisfy the additional cubic condition w = 0. Keeping u arbitrary, the only way 
to satisfy this constraint is to set 7=0; indeed, the inequalities (2.6) immediately imply: 
b = c = [3 = and a = a and therefore a vanishing w. The matrix 7i in (2.8) takes now a 
very simple form, and the non- vanishing eigenvalues are complex: = a±iu. Since 

a is positive, most of the entries of the evolution matrix M(t) in (3.2) are still exponentially 
suppressed for large t; however, the presence of the zero eigenvalue now implies M 33 (t) = 1, 
so that the asymptotic form of (2.11) changes: 

Vu.^{t) ~ ^ sin 2 20. (3.6) 

The large-time behaviors (3.5) and (3.6), for the particular case 7=0, are char- 
acteristic of the presence of the dissipative contribution (2.5) to the evolution equation 
(2.2). However, in general, it might be very difficult to distinguish these behaviours from 
the one obtained in the standard case. Although in principle P^L^ (£) in (2.13) has a 
purely oscillatory form, in any actual observational condition, the oscillations are likely 
to be averaged away, so that also in this case (3.6) holds. Therefore, when the mixing is 
maximal (sin 2 29 ~ 1), or in the special situation in which only one dissipative parameter 
is non- vanishing (7 = 0), the asymptotic large t behaviors (3.5) and (3.6) turn out to be 
indistinguishable from that of vi® L>v (t) . In these cases, one has to study the full time 
dependence of the transition probability. 



4. TRANSITION PROBABILITY: EXPLICIT FORM 

The general expression of the transition probability V Va ^ v it) in terms of u> and 
the dissipative parameters is very complicated and not particularly useful in practical 
applications. Therefore, we shall now discuss some approximations for which V Va ^ v (t) 
assumes a more manageable form; it might be of interest to compare these expressions 
with actual experimental data in order to put limits on the magnitude of the dissipative 
constants. Although this is clearly beyond the scope of the present investigation, we shall 
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nevertheless briefly comment about the rough sensitivity that one might expect from the 
analysis of present and future experiments. 

As discussed before, in general V V( ,^ V At) contains two kind of contributions: oscillat- 
ing terms, controlled by u>, and exponentially dumping terms, signaling dissipative effects. 
The relative dominance of these two types of behaviour depends on the magnitude of a, 6, 
c, a, (3 and 7 when compared to u>. 

In our approach, the dissipative contribution (2.5) to the evolution equation (2.2) 
should be regarded as phenomenological; it is therefore hard to give an apriori estimate of 
the magnitude of the non-standard parameters in L[p\. As mentioned in the Introduction 
and further discussed in the next section, a general framework in which dissipative effects 
naturally emerge is provided by the study of open quantum systems, i.e. systems in weak 
interactions with large environments. In such cases the dissipative effects can be roughly 
estimated to be proportional to the typical energy scale of the system, while suppressed 
by inverse powers of the characteristic energy scale of the environment. [1-3, 16, 7] 

In the case of the neutrino system, on the basis of these considerations and in line 
with the idea that dissipation is induced by quantum effects at Planck's scale, one expects 
the values of the parameters a, 6, c, a, (3 and 7 in (2.5) to be very small; for any fixed 
neutrino source and observational conditions, an upper bound on the magnitude of these 
parameters can be roughly evaluated to be of order E 2 /Mp, with Mp the Planck mass. 
The ratio of a, 6, c, a, (3 and 7 with u can thus be estimated to be at most of order 
10~ 10 E 3 /Am 2 , with E expressed in MeV and the neutrino square mass difference Am 2 
in eV 2 . By taking for E and Am 2 values that are typical of various neutrino sources, this 
ratio turns out to be about 10 2 for atmospheric neutrinos, of order one for solar neutrinos, 
while for accelerator neutrinos it can be as small as 10~ 2 . 

When the dissipative, non-standard parameters are large or of the same order of 
magnitude of u>, all entries of H in (2.8) are in general different from zero. In this case 
a useful approximation is to assume c and (3 to be much smaller than the remaining 
constants. ^ To lowest order, the matrix Ti becomes block diagonal and a manageable 
expression for the transition probability in (2.11) can be obtained. Explicitly, one finds: 



V Ue ^ (t) = l - 1 cos 2 26 [1 - e" 27t ] + sin 2 26 



1 - e~ At ( cos(2 Q t) + — e — sin(2 Q t) 

2 1 In 



(4.1) 



where 



A = a + a , B = a-a + 2ib , O = V^ 2 ~ l s l 2 / 4 • ( 4 - 2 ) 



The oscillating behavior in (4.1) depends on the magnitude of lo with respect to \B\; 
when lo < |S|/4, the frequency Qq becomes purely imaginary and V Ue ^ u ^(t) contains only 
exponential terms. In any case, the exponential dumping terms in (4.1) dominate for large 
t, and the limit (3.5) is recovered. 

A further simplification occurs when 7 = ; as already observed in the previous 
section, this automatically guarantees c = (3 = 0, an further imposes 6=0 and a = a. In 



t Note that this choice is perfectly compatible with the constraints of complete positivity 
given in (2.6) 
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this case, (4.1) reduces to: 



V„ e ^(t) = ±sm 2 20 



-2at 



cos(2 uji) 



(4.3) 



This is the most simple form that the transition probability formula takes in presence of 
dissipative effects: with respect to the standard expression in (2.13), (4.3) contains an 
exponential dumping factor in front of the oscillating term. It can be used to derive the 
rough order of magnitude bound on the non-standard parameter a that can be expected 
from neutrino experiments. Assuming that the dumping due to the exponential term is 
not exceeding a few percent, from (4.3) one derives: at < 1. Since the neutrinos are 
relativistic, the flight-time between emission and detection is roughly the same as the dis- 
tance £ between source and detector. Then, one has: a < l/£, where l/£ is approximately 
10~ 22 GeV, 10~ 24 GeV, 10~ 27 GeV for accelerator, atmospheric, solar neutrinos, respec- 
tively. Although the best bound on a seems to be given by solar neutrinos experiments, 
due to the larger £, atmospheric neutrinos data are the most suitable for a meaningful fit 
of (4.3), since in this case its time (or £) dependence can actually be probed. 

Another very useful approximation of the general formula (2.11) for the transition 
probability can be obtained when the non-standard parameters a, 6, c, a, (3 and 7 are 
small compared with uj. In this case, the additional dissipative term L[p] in (2.2) can be 
treated as a perturbation. Then, up to second order in the small parameters, one explicitly 
gets: 



V,, .„„(/) cos 2 2^ 



1 ( 1 si.r(0/) 



+ sin 2 29 
+ sin4# e~ At 



^ 2 

1 , w ( cos(2fit) + ^ s in(2fit) - sin 2 (fit) 

2 UD U u 



TZeC . ^, TZe[C(A + B - 2 7 )1 ,' 
sin(2 fit) + 1 V Q2 ^ sin 2 (fit) 



where A and B are as in (4.2), while: 



C = c + i/3 , fi = v 7 ^ 2 - |C| 2 - |S| 2 /4 



(4.4) 
(4.5) 



In the previous formula, we have reconstructed the exponential factors by consistently 
putting together the terms linear and quadratic in t; a similar treatment has allowed 
writing the oscillatory contributions in terms of the frequency flJ 

As a further check, note that the expression (4.4) reduces to that in (4.2) for \C\ = 0, 
i.e. when c = f3 = 0: it is therefore a correction to (4.2) for nonvanishing C. In this 
respect, the validity of (4.4) goes beyond the approximation in which it has been derived, 
since it can be considered as the expansion up to second order of the general formula (2.11) 
for c and (3 small. Therefore, it can be used with confidence in fitting experimental data 
from neutrino oscillation experiments. 



t This frequency is now real, since by hypothesis uj 2 3> \C\ 2 + |S| 2 /4. 
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In this respect, the data on atmospheric neutrinos are presently the best place to look 
for dissipative effects. Applying techniques similar to the ones employed e.g. in [28] and 
[29] to the generalized transition probability (4.4), one should be able to extract from the 
actual data useful bounds on some of the non-standard parameters in (2.5). Nevertheless, 
one should note that having in general six additional unknowns to fit will certainly make 
the procedure much more difficult and complex than in the standard case, where only the 
mixing angle 9 and the mass difference Am 2 are present; only for the simplified expression 
(4.3), that contains just one additional parameter besides 9 and Am 2 , one can actually 
expect a good fitting accuracy. 



5. DISCUSSION 

In the previous sections we have discussed how a phenomenological approach based 
on quantum dynamical semigroup can be used to describe dissipative dynamics for the 
neutrino system. As already mentioned in the introductory remarks, this phenomenological 
treatment can be supported by physical considerations. Indeed, a general picture in which 
dissipative effects naturally emerge is provided by systems in weak interaction with suitable 
environments. In the case of elementary particle systems, these effects are likely to originate 
from the dynamics of strings; however, an effective description of the environment, encoding 
some of the "collective" properties of the underlying fundamental theory, is quite adequate 
for a more physical discussion of evolutions of type (2.2), (2. 3). [7] 

To be more specific, in the case of neutrino systems, the total hamiltonian can always 
be decomposed as: 

H tot = H eff <g> 1 + 1 <g> H £ + gH' , (5.1) 

where H e Q is as in (2.4), while He describes the internal dynamics of the environment £. 
The interaction terms between the two systems are assumed to be weak: they are encoded 
in H' ', with g a small coupling constant. 

Furthermore, the mechanism of neutrino production is different from the one respon- 
sible for the dissipative effects; it is therefore natural to assume that the neutrino state 
and that of the environment be uncorrelated at the moment of the neutrino emission. In 
other words, the initial state of the total system can be taken to be in factorized form: 
(hot = P® p£- 

The time evolution of the neutrino state p, obtained by tracing over the environment 
degrees of freedom, 

p h-> p(t) = Tr £ L' 111 ^ 1 (p <g> ps) e lHtott ] , (5.2) 

is in general very complicated and can not be described explicitly. Nevertheless, an evolu- 
tion equation of the form (2.2), (2.3) for p(t) naturally emerges by taking into account the 
physical requirement that the interaction between neutrinos and environment be weak. 
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There are essentially two different ways of implementing this condition: [1-3] they cor- 
respond to the two ways of making the ratio t/t£ large, where r is the typical variation 
time of p(t), while T£ represents the typical decay time of the correlations in the environ- 
ment. Indeed, only for r ^> one expects the memory effects implicitly encoded in (5.2) 
to be negligible, and a local in time evolution for p(t) to be valid. 

When T£ becomes small, while r remains finite, one speaks of "singular coupling 
limit", since the typical time-correlations of the environment approach a 5-function. In 
the other case, when T£ remains finite, while r becomes large, one works in the framework 
of the so-called "weak coupling limit"; in practice, this is obtained by rescaling the time 
variable, t — > t/g 2 , and sending the coupling constant g to zero (van Hove limit). 

The choice between the two limits is made on the basis of physical considerations. In 
the case of unstable systems for instance, the weak coupling choice is unviable, since in this 
case r can be identified with the (finite) lifetime. On the contrary, for neutrino systems 
both limits are in principle allowed. ^ They give rise to different explicit expressions for the 
additional contribution L[p] in (2.3); in the case of the singular coupling limit, one finds: 

L[p] = -J o °° dt Tr e {[e iHet H' e~ iHet , [lf',p<g> p £ ]] } , (5.3) 

while in the weak coupling limit, one obtains: 

L[p] = -\un^^ J T ds J™ dt Tr £ ^e lH ^ s [e iHot H' e~ iHot , [H\ p <g> p s ]] e~ lH ^ s J , 

(5.4) 

where H Q is the limit of H tot when the coupling constant g vanishes. 

As mentioned before, the general form of the expressions for L[p] given above does 
not actually depend very much on the details of the environment dynamics; an effective 
description that takes into account its most fundamental characteristic properties is enough 
to allow an explicit evaluation of the integrals in (5.3) and (5.4). Following the idea that 
the dissipative effects originate from the low energy string dynamics at Planck's scale, one 
can effectively model the environment as a gas of DO-branes, in thermodynamic equilibrium 
at Planck's temperature; these quanta obey an infinite statistics. [30-32] 

Explicit computations then show that both expressions (5.3) and (5.4) assumes pre- 
cisely the form given in (2.5).$ However, while in the case (5.3) all six parameters a, 6, 
c, a, (3, 7 are in general nonvanishing, in the weak coupling limit the average procedure 
in (5.4) implies a = a and b = c = (3 = ^= 0, independently from the value of ui. As a 
consequence, when the weak coupling limit conditions are satisfied, the dissipative piece of 



t Nevertheless, it should be stressed that the condition that makes the characteristic 
times of the neutrino system much larger than that of the environment, implicit in the 
weak coupling limit, might not be attainable in all situations; on the contrary, the condition 
on the environment time-correlations necessary for the singular coupling limit seems more 
natural, in view of its possible "stringy" origin. 

* The steps followed for the evaluation of the integrals in (5.3) and (5.4) do not much 
differ from the ones presented in [7]; the details are therefore omitted. 
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the extended dynamics is controlled by a single parameter and the transition probability 
V Ve ^ v (t) assumes the simplified form presented in (4.3); on the other hand, the more 
general behaviour (4.4) is surely the result of a singular coupling limit procedure. 

Therefore, the indication of a non- vanishing value for more than one of the parameters 
a, 6, c, a, (3, 7 in neutrino oscillation experiments would certainly select the form (5.3) 
for the dissipative piece L[p]; in turn, this would provide some indirect informations on 
the structure of the environment and thus on the effective dynamics of low energy string 
theory. 

In closing, we would like to make a few comments on the existing literature on the 
subject. Studies of possible phenomena violating quantum mechanics in neutrino dynam- 
ics have recently appeared. [33-35] Based on ideas originally presented in [16], they discuss 
modifications of the standard oscillation probability formula. However, the extended dy- 
namics used in such investigations is that of [16], which does not satisfy the condition of 
complete positivity; as mention before, this could lead to serious inconsistencies. We stress 
that to avoid these problems, one has to adopt phenomenological descriptions based on 
the equations (2.2) and (2.5). 

Kinetic evolution equations similar to the one presented in Section 2 have been used 
to describe other, more conventional dissipative phenomena that arise due to the scat- 
tering and absorption processes in the core of supernovae or in the early universe. [36] In 
these extreme conditions, the frequent collisions affect the free evolutions of the neutrino 
species, and the consequent decoherence effects modify the oscillation pattern. The phys- 
ical situation is now quite different from the one discussed in the previous sections and 
necessarily requires a second-quantized, field-theoretical extension of the formalism. Fur- 
ther, the derivation of the evolution equations can not rely on the weak-coupling limit 
arguments discussed above; rather, it is based on the use of specific effective interaction 
hamiltonians. Nevertheless, also in these cases physical requirements like the condition of 
complete positivity should in general be enforced and might turn out to be crucial for the 
self-consistency of the formalism. 

Dynamical equations of the form (2.2) have further been employed for the study of 
the propagation of neutrinos in a density fluctuating media, in particular, in the interior 
of the sun. [37] They give rise to expressions for the surviving probability of the electron 
neutrinos that differ from those obtained in the framework of standard matter oscillations. 
Although described in terms of quantum dynamical semigroups, these density fluctuation 
have their origin in the dynamics of the sun and operate at energy scales quite different 
from Planck's mass. Therefore, they can be easily isolated from the dissipative effects 
discussed in the previous sections, that, in view of their "microscopic" origin, are not 
expected to be influenced by long-range phenomena. 

The recent experimental data, in particular on solar and atmospheric neutrinos, show 
evidence of attenuation in the expected neutrino flux, signaling disappearance phenomena. 
Although one is led to interpret these results in terms of the standard oscillation formula 
(2.13), several other physical mechanisms have been proposed as alternative explanation 
for the effect, in particular: neutrino decay, flavour changing neutral currents, violation 
of Lorentz invariance or of the equivalence principle. In all these cases, the transition 
probability V has a dependence on time (or pathlength) and neutrino energy that differ 
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from the standard one. (For recent discussions, see [28, 29].) 

The dissipative effects studied here are clearly distinct and independent from all these 
explanations for the neutrino flux deficit. In particular, the dependence of V on the non- 
standard parameters a, b, c, a, (3 and 7 is distinctive of dissipative phenomena and can 
not be mimicked by the other mechanisms. This is a great advantage in the process of 
fitting and comparing the experimental data, since it makes possible the identification of 
the dissipative contributions quite independently from all other effects. 



Note Added 

After the submission of our manuscript, the paper in Ref.[38] appeared; using the at- 
mospheric neutrino data of the Super-Kamiokande experiment a bound on one of the 
dissipative parameters was obtained. 
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